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We present the existence of the Kosterlitz-Thouless (KT) transition for n-atic tangent-plane 
order on a deformable spherical surface and investigate the development of quasi-long range n-atic 
order and the continuous shape changes of a deformable spherical surface below the KT transition 
in the low temperature (mean field) limit. The n-atic order parameter ip = exp[inO] describes, 
respectively, vector, nematic, and hexatic order for n = 1,2, and 6. We derive a Coulomb gas 
Hamiltonian to describe it. We then convert it into the sine-Gordon Hamiltonian and find a linear 
coupling between a scalar field and the Gaussian curvature on the fluctuating spherical surface. After 
integrating over the shape fluctuations, we obtain the massive sine-Gordon Hamiltonian, where the 
interaction between vortices is screened. We find, for n 2 K n /n 3> 1/4 where K n and k are n-atic 
j and bending rigidity, respectively, the KT transition is supressed altogether. On the other hand, 

if n 2 K n /n <C 1/4, there is an effective KT transition. In the ordered phase, tangent-plane n-atic 
order expels the Gaussian curvature. In addition, the total vorticity of orientational order on a 
surface of genus zero is two. Thus, the ordered phase of an n-atic on such a surface will have 2n 
vortices of strength 1/n, 2n zeros in its order parameter, and a nonspherical equilibrium shape. 
' Our calculations are based on a phenomenological model with a gauge-like coupling between tp and 

£ — , curvature, and our analysis follows closely the Abrikosov treatment of a type II superconductor just 

£S| ' below H C 2. 
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I. INTRODUCTION 

00 

When surfactant molecules, or amphiphiles which consist of molecules that combine both polar and non-polar 
^\ \ parts, are dissolved in a single solvent such as water, these molecules tend to form bilayer membranes where the 
"p-; • hydrocarbon parts of each monolayer are aggregated in the middle of the bilayer to reduce the contact between the 
water and the non-polar parts of the molecule due to the hydrophobic interactions of the hydrocarbon chains. These 
membranes are not covalenly bonded, but rather stabilized by the weaker hydrophobic interaction. These membranes 
can have characteristic sizes which may be much larger than those of a single molecule. The sizes and shapes of these 
membranes can change as a function of temperature, salinity, and/or surfactant concentration. Depending on the 
physical conditions of the system, these membranes can form random extended surfaces, regular periodic structures, 
[ or closed vesicles separating an interior from an exterior jlj. Furthermore, some bilayer membranes are prototypes of 
• • ■ biological systems, although it should be noted that true biological membranes have several components and greater 
complexity. 

In a wide class of membranes, molecules move freely within the membrane forming a 2-dimensional fluid offering 
little resistance to changes in membrane shape. Such membranes are physical examples of random surfaces, which 
can undergo violent shape changes. Molecules in membranes can also exhibit varying degrees of orientational and 
positional order including tilt(smectic-C), hexatic, and crystalline orderings similar to those found in free standing 
liquid crystal films These ordered membranes provide fascinating laboratories for the study of the coupling 

between order and geometry analogous in many ways to the coupling between matter and geometry in general 
relativity. The underlying cause of this coupling is as follows. A field describing orientational order cannot be parallel 
everywhere if it is forced to lie on a surface, such as that of a sphere, that curves in two directions and thus has 
nonzero Gaussian curvature. It could lower its energy by flattening the surface. An additional complication arises 
when order develops on closed surfaces. A closed surface can be classified according to its genus g: the number of 
handles. Orientational order on a closed surface necessarily has topological defects (vortices) with total strength 
(vorticity) equal to the Euler characteristic 2(1 — g) of the surface ||. Tangent-plane order on a sphere (torus) will 
have vorticity 2 (0), since a sphere (torus) has genus (1), respectively. The continuous development of vector order 
on a deformable surface of genus zero will be accompanied by a continuous change from spherical to ellipsoidal shape 
||. Since vortices are energetically costly, it may be favorable for a closed physical membrane to transform into an 
open cylindrical structure when tangent-plane order develops in response to changes in temperature or other control 
variables . Indeed, there are a number of experimental examples of shape changes that may be explained by the 
development of tangent-plane order. 
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Orientational orders for smectic-C and hexatic liquid crystals are described by an n-atic order parameter ip — 
(exp[in6]) with n = 1 and n — 6, respectively. In this paper, we present the existence of the KT transition of n-atic 
order on a closed surface of genus zero and investigate the development of n-atic order below the KT transition 
temperature and the concomitant change in shape from spherical to non-spherical. An n-atic order parameter can 
have vortices of strength l/n, and, since it is generally favorable to form vortices of minimum strength, we expect 
2n maximally separated vortices of strength l/n to be present in the ordered phase (low temperature phase). Thus 
for n = 1, we expect two antipodal vortices, and respectively for n = 2, n = 3 and n — 6, we expect vortices at 
the vertices of a tetrahedron, an octahedron, and an icosahedron. Indeed calculations on a rigid sphere confirm this 
conjecture ||. We also expect the vesicle shape to change from spherical to ellipsoidal, tetrahedral, octahedral, or 
icosahedral in the four cases above. 

Our calculations are based on a phenomcnological Hamiltonian for a complex order parameter field whose coupling 
to shape occurs via a covariant derivative and via changes in the metric tensor ||. The model is almost identical 
to the Landau-Ginzburg theory of superconductivity except that vorticity is fixed by surface topology rather than 
energetically determined by an external magnetic field. The ordering transition we find is very similar to the transition 
from a normal metal to the Abrikosov vortex lattice in a superconductor, and indeed our analysis follows very closely 
that of Abrikosov JToj ] . We find the order parameter, which is necessarily spatially inhomogeneous because of the 
topological constraint on the total vorticity, from among a highly degenerate set of functions that diagonalize the 
harmonic Hamiltonian on a rigid sphere. This degenerate set has exactly 2n zeros at arbitrary positions on the sphere 
and is very similar in form to fractional quantum Hall wave functions 

This paper is organized as follows. In Sec. II, we review fluid membranes in some detail and describe their free 
energy as a function of curvature. We introduce a tangent-plane orientational order parameter and its coupling to 
shape changes in Sec. III. The model Hamiltonian for a deformable spherical membrane with n-atic order is presented 
in the normal gauge representation and the existence of the KT transition is described in Sec. IV. Minimization 
scheme and the calculated shapes for temperatures below the KT transition temperature to the ordered n-atic phase 
for n = 1, 2, 3, 4, and 6 are given in Sec. V, and we give concluding remarks in Sec. VI. 

II. FLUID MEMBRANES 

In this paper, we use the term membrane to denote a thin film of one material that separates two similar mate- 
rials (bilayer membrane). We focus on fluid membranes and their possible degrees of orientational order. In fluid 
memebranes, there is no in-plane shear modulus and the only in-plane deformations are compressions/expansions. 
An n-atic membrane, although fluid, will sustain a local orientational order defined modulo rotations by 2n/n, which 
is described by a local order parameter ip = exp[m0]. 

Although fluid membranes can be composed of many different types of chemical and molecular species, their 
behavior (shapes, fluctuations, thermodynamics) can be understood from a unified point of view that considers their 
free energy of deformation. If the membrane were constrained to lie in a plane, the only relevant energy would be 
the compression of the molecules, that is, change of the average area per molecule. This is analogous to sound waves 
in a 3-dimensional fluid; there is no low-frequency reponse of the system to shear. However, since the membrane 
can also deform in the normal direction, there is an additional set of modes describing the conformations of the film. 
These out-of-plane deformations are known as bending or curvature modes and the free energy associated with such 
modes is known as the curvature free energy. While a general deformation of the membrane involves changes in 
both volume and curvature, we shall see that the lowest energy deformations usually involve only curvature. In most 
systems, changing the average volume of the membrane is a higher energy process and hence is less important when 
effects involving the thermal behavior of the membrane are considered. For a membrane with finite thickness, we 
denote as pure curvature deformations those perturbations of the membrane that do not change the overall membrane 
volume, but where there may be local stretching and compression of different parts of the film. The location of the 
interface within a membrane can be chosen to lowest order in the curvature energy so that the surface which defines 
the interface undergoes no stretching or compression (neutral surface). For long- wavelength curvature deformations 
whose wave lengths are much larger than the membrane thickness, the exact position of the neutral surface within 
the membrane is not crucial. 

The membrane is characterized by the area per molecule, E, its thickness, A, and its curvature. For simplicity, we 
assume that the equation of state of the membrane determines the thickness as a function of the area per molecule. 
A simple example is the case of an incompressible molecule where the product of AS is constrained to equal the 
molecular volume so that A ~ I/E. We thus take the curved membrane to be characterized only by the area per 
molecule and its curvature. Consider a locally flat, isolated interface. Saturation occurs when the free energy achieves 
a minimum with respect to E; dfo/dY< = where /o is the free energy per molecule for a flat layer and E is the area 
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per molecule. The free energy per molecule is minimized when E = Eo, the optimal value of the area per molecule 
which arises from a balance of terms such as the entropy and the tension terms or attractions. The entropy favors a 
larger area per molecule because of the greater number of center-of-mass positions and chain conformations, while the 
tension terms and attractions favor a small value of E to reduce the contact of the hydrocarbon chains with the water. 
There can be deviations in the area per molecule from this optimal value and the energy cost of such a compression 
or expansion is 

A/o^/^E-Eo) 2 (2.1) 

where the primes signify a derivative with respect to E. However, these deformations are typically of higher energy 
than the curvature deformations. A membrane can change its shape or size with a much lower free energy cost than 
that required to compress or expand it. It is important to remember, therefore, that for insoluble amphiphilcs it is 
the saturation of the membrane and the minimization of the area per molecule that permits the usual surface-tension 
term to be neglected; df/dT, = 0. The surface tension is no longer relevant since the molecules adjust their area to 
optimize the free energy and it is the curvature energy that determines the properties of the membrane. 

To construct the effective free energy for fluid membranes, it is convenient to introduce general curvilinear coor- 
dinates u = {w 1 ,^ 2 } and the metric structure. Defining locally a system of coordinates u on the membrane and 
denoting by X(u) the position of the point u in bulk <i-dimensional Euclidean space, the metric tensor g a p(u) induced 
by the embedding is 

<Mu) = t„(u) • t^u) ; t a = (2.2) 

where t Q is a tangent vector to the surface, and the element of area is 

d 2 s = d 2 u^fg ; g = det(g af3 ). (2.3) 

The extrinsic curvature tensor K Q ^ is defined by 

K Q/3 (u) = D a D X(u), (2.4) 

where D a is a covariant derivative with respects to the metric g a p- In general, K Q £j is normal to the surface [3]. 
Therefore, in the particular case of a surface in i? 3 , K Q ^ is proportional to the unit normal vector N to the surface 
and it is written as 



K a/3 = K ap -N, (2.5) 

where K a p is a symmetric tensor. We can dicuss the curvature energy using symmetry considerations. The free energy 
must be only a function of the field X(u), invariant under displacements and rotations in R d , and reparametrization- 
invariant. Expanding in local terms involving more and more derivatives and keeping only the terms relevant by naive 
power counting, the most general form of the curvature energy TL C up to quadratic order in the curvatures can be 
written in terms of the mean and Gaussian curvatures and has only two terms 

n c = ^nj d 2 u^(H - H ) 2 + hj d 2 u^K. (2.6) 

Here, H is the mean curvature, \K^ : and K is the Gaussian curvature, detiff . This form for the curvature energy 
was discussed by Helfrich and Canham and states that the mean curvature that minimizes the energy has a value 
Hq, termed the spontaneous curvature of the membrane p2| . The energy cost of deviating from the spontaneous 
curvature is the bending or curvature modulus, k. The parameter R 7 known as the saddle-splay modulus, measures 
the energy cost of saddle-like deformation. The spontaneous curvature describes the tendency of the bilayer membrane 
to bend. It is viewed to arise from the fact that the two layers of the bilayer do not have areas per molecule which 
are exactly equal, nor are their curvatures equal and opposite. This can result in a vesicle that is more stable than a 
flat, lamellar bilayer since it may be that the favorably curved outer layer has a smaller value of the area per molecule 
and hence more molecules in it than the inner layer. The bending moduli, k and R, arise from the elastic constants 
determined by the head-head and tail-tail interactions. It is expected that these moduli are senstively dependent on 
the surfactant chain length but only weakly dependent on the head-head interaction strength. The parameters Hq, k, 
and R can be derived from a simple microscopic model that incorporates both the change in the area per molecule and 
the curvature. We note that a stable membrane will always have k > 0. However, the sign of R can be either positive 
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or negative. Membranes that prefer isotropic shapes where the Gaussian curvature K > such as spheres or planes 
will have re > 0, while membranes that prefer saddle shapes where the Gaussian curvature K < will have re < 0. 
One can show that the requirement of a positive-definite quadratic term implies that membranes are only stable if 
re — 2re > 0; otherwise higher order curvature terms are needed to stabilize the system. We consider the spherical 
membrane whose topology is fixed. Thus the Gaussian curvature energy gives constant contribution by Gauss-Bonnet 
theorem J ^fgK — 2irx- 

III. ORIENTATIONAL ORDER 

In a fluid membrane, molecules can flow freely to adapt themselves to any particular shape of the surface. If 
correlations among the molecular positions of the molecules forming the membrane exist, the molecules may exhibit 
in-plane crystalline order and form a kind of two-dimensional solid. On the other hand, the molecules may exhibit 
a weaker order in which the orientations are correlated at long distance scale. Orientational order means that to 
each point on the membrane is associated a preferred direction whithin the tangent plane of the membrane. For 
example at high temperature, the stable phase of the membrane is generally the liquid phase with no translational or 
orientational order if we are considering the hexatic order, or the smectic-A phase in which molecular axes are normal 
to the surface for the vector order. At lower temperature, the membrane can condense into an hexatic phase in which 
there is quasi-long-range 6-fold bond-angle order or into a smectic-C phase in which molecules tilt relative to the 
surface normal. To describe Sm-C (vector) and/or hexatic order, we introduce at each point X(u) on the membrane 
a unit vector m(u) in the tangent plane of the membrane. For Sm-C order, m(u) is a true vector, invariant under 
rotations of 2pir (p is an integer) about the unit surface normal N(u) erected at X(u). For hexatic order, rotations 
of m(u) by 2pir/6 about N(u) lead to physically equivalent states. More generally, we consider "n-atic" order in 
which rotations of m(u) by 2pn/n produce physically equivalent states. A two-dimensional nematic with an in-plane 
symmetric traceless tensor order parameter is an example of a 2-atic. Although we know of no physical realizations of 
other n-atics yet, we find it instructive to consider how the development of such order affects morphological changes 
in spherical vesicles. 

To describe tangent-plane n-atic order, we introduce orthonormal unit vectors ei and e-i at each point on the 
membrane. ei(u) • m(u) = cos 0(u) defines a local angle 0(u). n-atic order is then described by the order parameter 
ip(u) = exp[inO(u)], which can be related to the n-th rank symmetric traceless tensor constructed from the unit vector 
m. The n-th rank symmetric traceless tensors are the n-th rank spherical tensors. They are listed in Table [| for 
n = 1, 2, 3, 4, and 6. In general, there are 2™ components in the n-th rank spherical tensor Q$ $ ...j since i±, 12, • • • , i n 
can be either 1 or 2. By permutational symmetry, there are only (n+ 1) possible independent components. However, 

there are (n — 1) additional traceless conditions. Hence, there are only 2 independent components in ^— « ■ The 
linear combination of these n-th rank spherical tensors 

E *' fe Q?l~i„' (3-7) 

il,i2,---,in 

where k is the number of 2's in (ii, 12, ■ ■ ■ , i n ), becomes (mi+imj)" where mi = ei-m = cos and 1712 — e2-m = sin0. 
Thus n-atic order parameter is described by 

i> = E i fe Q£L, n = e ine (3-8) 

Note that since 0(u) depends on the choice of orthonormal vectors ei and e2, the order parameter tf>(u) does as well. 
This means that any spatial derivatives in a phenomenological Hamiltonian for ijj must be covariant derivatives. 

We now describe the Hamiltonian of the membrane with orientational order in a reparametrization-invariant way. 
For a flat membrane the vector Hamiltonian corresponds to the usual XY model 

H = T^Kn J d 2 xd a m ■ d a m. (3.9) 

Using the fact that any spatial derivatives must be covariant derivatives and 

D a m = (D a m )t p (3.10) 
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is the tangential component of d a m, one can show that the only possible free energy term for m which respects this 
Z\ symmetry is 



H = ^Ki J d 2 u y /g'D tl m a D f *m a . (3.11) 



The dimensionless coupling constant K\ is the vector-order stiffness and measures the strength of the coupling between 
the orientations of neighboring vectors. We then generalize this XY-like Hamiltonian for the vector order to the n-atic 
Hamiltonian using the n-th rank spherical tensors introduced in Table [|. In general, the n-atic Hamiltonian can be 
written as 

U = ±K n / dVff ^ai „ ai "' an , (3-12) 

2 J V Q(«) 1 n Q(«) Q a 

where K n is the n-atic rigidity. Using the spherical tensors in Table 1, we find 

4 4 ttl " 2(n-l)' 

and 

2 

^QW^-^^QW^...^ = -JL— Dfl m a D>*m a , (3.14) 

for n > 1. Hence, the Hamiltonian for the n-atic order becomes 

H n = \n 2 K n J d 2 u^D^m a D»m a . (3.15) 

In this form of TL n , we neglected all terms that are irrelevant at large distance by power counting. Other terms 
such as m a KoKP' y my which couple m to the principal directions of curvature of the membrane are not invariant 
under the global rotation by 27r/n of m, and are irrelevant at large distance. Thus, this free energy has a full 0(2) 
rotational symmetry. This is similar to the fact a 2-dimcnsional crystal with hexagonal or triangular structure has 
isotropic elastic properties at large distance scale. For n-atics with n > 3, there is only one elastic constant K n . For 
n = 1 or n = 2, there are in general two elastic constants. For simplicity, we will consider the single elastic constant 
approximation for all n-atics. 

To describe this Hamiltonian in terms of a local angular order parameter 0, at each point we introduce two or- 
thonormal vectors e tt (u)(a = 1, 2) tangent to the membrane. This is equivalent to introduce a 2-bein e l a (u) compatible 
with the induced metric g a p (u). In components, 

e a = e a a t a (3.16) 

and the orthonormality e a • = 8 a b implies 

<%49*fi = Sab; e^e b 5 ab = g^. (3.17) 

The angular order parameter is frustrated by the rotation of tangent vectors that occurs under parallel transport 
on a curved surface. The amount of frustration is given by the gauge field A a , i.e. the covariant derivative of e a 
in direction a defines the gauge field A a . Under parallel transport in direction du a , each e a is rotated by an angle 
A a du a . Thus the gauge field A a is defined by 

D a e a = -A a e ab e b , (3.18) 

where e a b is the antisymmetric tensor with £12 = —£21 = 1 and A a e a b is called the spin-connection and describes how 
the basis vector e a rotates under parallel transport according to the Gaussian curvature K of the surface. In fact, A a 
is related to K. The curl of the gauge field A a is the Gaussian curvature; 

y af3 D a A = K, (3.19) 

where 7"^ is the antisymmetric tensor defined via 
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(3.20) 



The covariant derivative of m = cos 9ei + sin 0e2 writes 



D a m = (D a m a )e a + m a (D a e a ) 
= (D a m a )e a - m a A a e ab e b 



= (D a Q)(— sin 9ei + cos 9e 2 ) 



— A Q (cos6e 2 — sin9ei) 
= (D a G- A a )ia x , 



(3.21) 



where m j 



sinOei + cos0e2 satisfying m • = 0. Then the n-atic Hamiltonian writes 



H n = -n 2 K n j d 2 u^gg al3 (d a Q - A a )(d f3 Q - A p ). 



(3.22) 



This form of the free energy is invariant under local transformations 0(u) — > 0(u) + A(u), A a (u) — ► A a (u) + d a A(u). 
This gauge invariance corresponds to a local rotation of the reference frame e a . In terms of the complex order 
parameter ip, the Hamiltonian for n-atic order becomes 



where D a i\) = (d a — inA a )?p. In this description, all n-atics have the same long-wavelength elastic energy. Their 
properties can differ however, because their topological excitations are characterized by different winding numbers. 
Since we assume that in the disordered state the membrane forms a spherical shape and there is no topological 
deformation of the membrane shape, order in -0 is necessarily accompanied by the topological excitations called 
vortices. As we can see in Ref. ||, the total winding vorticity of a vector field on a closed surface with h handles 
must be 2(1 — h). Thus, a vector field on the surface of a sphere (h = 0) has total vorticity 2. The minimum winding 
number disclination for an n-atic is 1/n. The energy of an individual disclination on both fiat and curved surfaces 
is proportional to the square of its winding number. It is therefore always favorable to form disclinations with the 
lowest possible winding number. In addition, disclinations with the same sign repel each other. These considerations 
imply that the ground state of a sphere with surface n-atic order will have 2n maximally separated disclinations of 
winding number 1/n. 



In this section, we present that in a certain region in the parameter space (re, K n ) there is an effective KT transition. 
Ovrut and Thomas have discussed the structure of the KT transition of a vortex-monopole Coulomb gas on a rigid 
sphere and have shown it is the same as in the planar case. They have shown the KT transition temperature on a 
rigid sphere is the same as that on the Euclidean plane, i.e. T^ eie — T^ nc — itKa/2, where Ka is the square of the 
vortex charge in their discription |l3| ]. Recently, we have extended the study of the structure of the KT transition 
to the hexatic Hamiltonian on the fluctuating spherical surfaces and have shown the hexatic Hamiltonian under goe s 
the effective KT transition at the same tenperature as in the planar case ]l^| . Following the procedure in Ref. [|l4[ , 
we generalize a study of the KT transition to the n-atic Hamiltonian and present the region in the parameter space 
(re, K n ) where the KT transition exists. 

From the previous sections, the n-atic Hamiltonian on the fluctuating surfaces writes TL = 7i c + H n ; 



We investigate the effect of thermal shape fluctuations of a genus zero surface on the KT transition in the low 
temperature limit (3k S> 1. In this limit, we can parametrize the surface by its radius vector as a function of standard 
polar coordinates u = (8, (j>) = f2; 




(3.23) 



IV. EXISTENCE OF THE KT TRANSITION 




(4.24) 



(4.25) 



G 



R(fi) = R(l + p(Sl))e r , 



(4.26) 



where e r is the radial unit vector and /o(f2) measures deviation from sphericity with radius R. This parametrization is 
a "normal gauge" . To make the equations simple, we map this parametrization onto the unit sphere parametrization 
with R = 1. Since we are interested in the limit f3n 3> 1 and in this limit Tl c dominates, we first minimize TL C over 
the shape fluctuation field p which gives p(Sl) = and then we minimize Tin over with p(Sl) = and find 



<59(fi) 



e=e° VS° 



X - dh g« a \d a Q« A° a ) = 0, (4.27) 



where the superscript stands for the rigid sphere with p(Sl) — 0. In Ref. ||, Lubensky and Prost show that in the 
ground state 2n disclinations of strength 27r/n are arranged at the vertices of polyhedron inscribed in the sphere. A 
disclination at u — with strength qi gives rise to a singular contribution 8q II1S to 0° satisfying 

du a d a ef le = q t , (4.28) 

r 

where T is a contour enclosing U;. Thus, in general d a Q° = d a <d' + where O' is nonsingular, u° = 9 Q 0Q lng , and 

^D a v° = 8° (SI), (4.29) 

where 

s°(n) = — y^6(Cl-Sk), (4-30) 

i=l 

where is the disclination density in the ground state and Sl^s are the coordinates of the vertices of icosahedron. Since 
(d a Q° — A° a ) satisfies D a (d a Q° — A°) = 0, it is divergence-less and purely transverse. Accordingly (d a Q° — A° a ) can 
be written in terms of the curl of scalar fields and by applying the operator ^ a Dp to (d a Q° — A Q a ), we find these 
scalar fields to be related to the Gaussian curvature K of the rigid sphere and the ground state disclination density 
on the rigid sphere, 

^ a D p (d a Q° Al) = j^D p v° a - j 0a D p A° a 

= s° - K , (4.31) 

where K is a Gaussian curvature of the rigid sphere and s° is the disclination density in the ground state. 
Now taking into account the bond angle fluctuations around 0° and the shape fluctuations around the sphere, 

= 0° + 0, A a = A° Q +SA a , (4.32) 

the full Hamiltonian writes Ti = Hq + 5H 

Ho = \n 2 K n J dSl(d a e Q - A° a )(d a e° - Al) 

m = \ K j dn { { y 2 + 2 )pf 

+ ^n 2 K n J dSl(d a Q - SA a ){d a e - SA a ) + 0(p 3 ). (4.33) 

The angle fluctuation field Q(Sl) can also have disclinations of strength q — 2n(k/n) where k is an integer, due to the 
thermal fluctuation Jl5| . Thus, d a Q can be decomposed into singular and nonsingular parts; <9 Q = <9 Q 0" +v a where 
©II is nonsingular, v a — d a Q smg and 

^DaVff = s{si), s(n) = 5^? i *(n-n i ), (4.34) 

i 

where s(Sl) is the thermally-excited disclination density with disclinations of strength qi at Sli. The vector v a can 
always be chosen so that it is purely transverse, i.e. D a v a = 0. In the hexatic Hamiltonian, d a Q always occurs in the 
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combination (d a Q — SA a ). The spin-connection SA a can and will in general have both a longitudinal and a transverse 
component. However, one can always redefine 0" to include the longitudinal part of SA a . This amounts to choosing 
locally rotated orthonormal vectors ei(u) and e2(u) so that D a SA a = 0. Thus we may take both v a and 5A a to be 
transverse and the hexatic Hamiltonian 



^n 2 K n J dn(d a e ]l + v a - SA a )(d a Q l] + v a - 5A a ) 
= Ti.\\ + 7i±, 

can be decomposed into a regular longitudinal part, 
and a transverse part, 

H x = \n 2 K n J drt(v a - SA a ){v a - 6A a ), 



(4.35) 



(4.36) 



(4.37) 



where the cross term J dfl(v a — 8A a )d a Q^ is dropped since D a (v a — A a ) = 0. 

It costs an energy e c (fc) to create the core of a disclination of strength k. (We assume for the moment that the core 
energies of the positive and negative disclinations are the same. See, however, Refs. fl^-|l8|.) Thus, partition sums 
should be weighted by a factor y^ = e~^ e "^ for each disclination of strength k. Since e c (fc) ~ fc 2 , we may at low 
temperature restrict our attention to configurations in which only configurations of strength ±1 appear. Let N± be 
the number of disclinations of strength ±1 and let u^± be the coordinate of the core of the disclination with strength 
±1 labeled by i. The hexatic membrane partition function can then be written as 



Z(K,K A ,y) = Tr v y N J VR J T>Qh- fm * e -K' H i +7i J, 



(4.38) 



where y = y%, and N = N + + N-. H± depends on all of the disclination coordinates Vl v ± where v — 1, 2, • • • , N±, 
and Tr v is the sum over all possible disclination distribution with the topological constraint |5J; 



Tr v = 



E 

N + ,N- 



nJnJ. 



n 



dn, 



n 



dn,. 



(4.39) 



where a 2 is a molecular solid angle. The Kronecker factor 8n + ,n_ in Tr v imposes the topological constraint that 
the total disclination strength on a sphere is 2 since with N + = we have 2n ground state disclinations with the 
strength 1/n giving the total discl ination strength 2n x (1/n) = 2. 

The hexatic model of Eq. (4.38) can easily be converted to a Coulomb gas model using 



-) aP D a {vp-SA p ) = s-5K = C, 



(4.40) 



which follows from Eq. (3.19) and Eq. (4.34) where 5K is the deviation of the Gaussian curvature from the rigid 
sphere. The quantity C = s — SK is a "charge" density with contributions arising both from disclinations and from 
Gaussian curvature. Equation (4.40) implies 



v a - SA a = -^^Dp—C, 



(4.41) 



where we used r y a \D x ^f a ^D a — — A and A = D a D a = (1/ \fg)d a ^fgg a P dp is the Laplacian on a surface with metric 
tensor g a p acting on a scalar. Recall [Eq. ( 3.20 )] that 7 Q ^ rot ates a vector b y 7r/ 2 so that (v a — 5A a ) is perpendicular 
to Dp(— A) _1 C and is thus manifestly transverse. Using Eq. ( f4.4l| ) in Eq. (1.37), we obtain 



Z = Tr v y N / VTL / VQUe-W.-flHi-fm. 



(4.42) 



where 
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He = \n 2 K n J dn^^c(n)^^c(n) 



(4.43) 



is the Coulomb Hamiltonian associated with the charge C. Since the longitudinal variable 8" appears only quadrati- 
cally in H\\ , the trace over 0" can be done directly giving the Liouville action |l9[ | arising from the conformal anomaly; 



VQh- 0n " = e~ fmL , (4.44) 



where 



f3H L = — U f dfl-— f dndtt'K(n) ( ~8(tl - fl')) K(fl'). (4.45) 
Sua 2, J 247r J \ A / 

The Coulomb gas partition function can thus be written 

Z = Tr v y N J me-^-^-^ . (4.46) 

The Coulomb gas model can be converted following standard procedures into a sine-Gordon model. The first step is 
to carry out a Hubbard-Stratonovich transformation on /37ic : 



where the Liouville factor e^ L is needed to ensure that e @ nc be one when C = 0. Inserting this in Eq. (4.46), we 
obtain 

Z = Tr v y N f PRMe-^-^e i / ,i!1(, - w)$ , (4.48) 

where 

/3W$ = ^(Pn^Kn)- 1 J dfld a <i>d a <i>. (4.49) 
The only dependence on disclinations is now in the term linear in n. Thus to carry out Tr v , we need only to evaluate 
Tr v2 /V/ df2s * 

N + / r \ N 

- S N m v n++n - I I — e 2ni *W n 1 

N + ,N. 



V - 5m m y N ++ N - ( [ — e ^^)/A + ( f ^ e -27ri*(n)/n" 

.fj>_N + lN-r N+ ' N - V U « 2 J \J a 2 

^ N+lN-l J 2ir \ y J a 2 

N + ,N- + v 

^ e (2y/a 2 ) / dnco S [27r(*/n)-w] ^ ^ 



Thus 



Z = J — J T>$ J VRe^ l3HK e^ l3n ' s, e {2y/a2) S <*ncos[27r(s/n)-w] e -i J dn$8K ^ 
We can now change variables, letting $ = (n/27r)(<j>' + w). The term linear in the Gaussian curvature then becomes 

- i [ dfl5K^-(cj + $') = { dn<$>'8K, (4.52) 

J 27T 27T J 
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where we used J dflSK — 0. The integral over uj in Eq. (4.51) is now trivial, and dropping the prime we obtain 



Z = / £>$ / 2?Re 



(4.53) 



where 



2y 



dfl cos $ — i- 



2tt 



(4.54) 



is the sine-Gordon action on a fluctuating surface of genus zero. The first two terms of this action are the gradient and 
cosine energies present on a rigid sphere. The final term provides the principal coupling between $ and fluctuations in 
the metric. It is analogous to the dilaton coupling po[| of string theory though here the coupling constant is imaginary 
rather than real. Note that the Liouville action is not explicitly present in Eq. (4.53). 

In the regime /3/t 3> 1, we can truncate the higher order terms in p. In the normal gauge, the partition function 
becomes 



Z = / VpV<$> exp 



~p K J dfi((V 2 + 2)p) : 

i/?r J d«(v$) 2 + ^| J dn cos$ 

i— I dfi$(V 2 + 2)p 
2n J 



(4.55) 



where (3T = l/(ATT 2 (3K n ) and we used 5K = (V 2 + 2)p. To lowest order in p, the shape fluctuation field p is linearly 
coupled to the scalar field $ which is the conjugate field to the disclinations. In Ref. |l(|, we have shown the similar 
coupling in the fluctuating flat membrane is quadratic in the shape fluctuation field. 

Integrating over the shape fluctuation field p gives the effective Hamiltonian for the conjugate field to the disclina- 
tions 



Z 



X>$exp 

~ I </<>C<>K<I> 



i/?r J dn ((v$) 2 + pH 2 ) 



(4.56) 



with p, 2 = n 2 K n / k. This is the massive sine-Gordon theory. The shape fluctuations induce the mass term for $ field 
and screen the Coulombic interaction between the disclinations giving the Yukawa interaction between them. This 
partition function is equivalent to that of the Yukawa gas Hamiltonian on the rigid sphere with radius R; 



H 



Yukawa 



(4.57) 



where 



G(n i -fi,-)=E 



21 + 1 



l(l + l) + p 2 



Pi (COS U)ij) 



P„{— COSCJy), 



(4.58) 



where P v (-) is the Legendre polynomial with degree v = — 1/2 ± (y/l — 4p 2 )/2 and Wjj is the angle between two 
disclinations at fii and ilj. For < p < 1/2, degree of Legendre polynomial, v, is real and the length scale introduced 
by \d = (p/R)" 1 is larger than the system size, 2R after recovering the original length scale by mapping the unit 
sphere back to the sphere with the radius R. On the other hand, if p > 1/2, v is a complex number, v = — 1/2 ± it 
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where r = {\J <±p 2 — l)/2 and Ad < 2_R. The length scale introduced by Ad — R/p may be interpreted as the Debye 
screening length arising from shape fluctuations. 

The interaction energy between two disclinations i and j at positions fli and Qj with strength and qj is given 
by qiqjG(dij) where dij = 2i?sin(o>jj/2) is the chordal distance between two disclinations on the sphere with radius 
R. The interaction G(dij) has the following limiting forms: 

C -i In {dij/2) , p' 1 > 2, rfy < 2i? 
G(dij) ~ < -| In {dij/2) , ii- 1 < 2, ^ < A d (4.59) 
[ e" d «/^, /i- 1 < 2, dy > A d 

Following the analogy of the 2D Coulomb gas, when the screening length is much larger than the system size, 
Xd 3> 2i?, the induced mass term arising from shape fluctuations is irrelavent for the KT transition and the KT 
transition temperature is given by Tq = nK n /2 for p <C 1/2. However, for p 1/2, the screening length is shorter 
than the system size, Ad <C 2R, and the mass term is relavent for the KT transition and changes the universality 
class of the system. There is no KT transition at non-zero temperature. The disclinations are always unbound at 
non-zero temperature and the KT transition temperature vanishes. Thus, when n 2 K n / n < 1/4, the effect of shape 
fluctuations is irrelavent and the KT transition occurs at the finite temperature. 



V. SHAPE CHANGES BELOW T c 

For n 2 K n /n <C 1/4, there exists a KT transition at the finite temperature Tc = nK n /2. To describe the phase 
transition and the corresponding shape changes, we introduce the magnitude of the complex n-atic order parameter 
as in ip — ipoe 1 and consider the simplest long-wavelength Landau- Ginzburg Hamiltonian for the n-atic order 
parameter ip; 

H LG = J d 2 u^ (V|^| 2 + \uW^j + \K n f d 2 u^g aP {d a -inA a )^{d [j +inAp)r- (5.60) 

This Landau-Ginzburg n-atic Hamiltonian is similar to the Landau-Ginzburg Hamiltonian for a superconductor in an 
external magnetic field, 



Hlg = / d 



73 „ 



r#| 2 + C|(V - ie*A)ip\ 2 + \u\^ |l(VxA-H) 2 

2 o7T 



(5.61) 



where e* = 2e/Kc. Both have a complex order parameter tp with covariant derivatives providing a coupling between 
ip and a "vector potential" A or A a . In a magnetic field, the superconductor can undergo a second order mean-field 
transition from a normal metal to the Abrikosov vortex lattice phase with a finite density of vortices determined 
energetically by temperature and the magnetic field H. The magnetic field is conjugate to the vortex number N v 
since J d 3 x(V x H) = LN v (f>o, where L is the length of the sample along H and </>o = hc/2e is the flux quantum. 
On a closed surface with n-atic order, there is a second-order mean-field transition to a state with vortex number 
determined by topology rather than conjugate external field. Thus, the mean-field n-atic transition on a closed surface 
is analogous to the transition to an Abrikosov phase with a fixed number of vortices rather than fixed field conjugate 
to vortex number. However, a Meissncr phase with zero vortices does not exist because we do not have an analogue 
for magnetic intensity H, rather topology fixes vortex number. 

The complete Hamiltonian for n-atic order on a deformable surface is H = H K + Hlg- The field p(Sl) in R(O) = 
i?o(l + p(^)) e r can be expanded in spherical harmonics. Any isotropic change in R can be described by Rq. In 
addition, uniform translations, which change neither the shape nor the energy of the vesicle, correspond to distortions 
in p with I = 1 and can be absorbed by fixing the center of mass of the vesicle. These considerations imply that p 
will have no I = or I = 1 components in spherical harmonic expansion: 

oo / 

p(«) = EE^™- ( 5 - 62 ) 

1=2 m=-l 

The shape and size of the vesicle are determined entirely by the parameters R and p\ m . The reduced tensors g a $ 
and K a j3 do not depend on Rq. Therefore, Rq can be expressed as a function of A and pi m via the relation 
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A = j d 2 Uyfg = R l j d 2 u^/g 



Hi j <i 



(l + P) 2 + 77 Wop) 



sin 2 9 



o(p 4 ) 
-o( P 4 ) 



(5.63) 



In the disordered phase above Tq , p — and R = Rq. We will use the Hamiltonian H. expressed in terms of reduced 
parameters and the constant area A in our calculations of shape changes below the second-order disordered-to-n-atic 
transition. 

Using the spherical polar coordinates with origin at the center of mass of the vesicle, the metric tensor g a p can be 
written as 



gap — d a R ■ d@TL — R-lga/3 
g a p = (l + P?g% + {d a p){dpp) 




9a P ~ \ sin 2 9 

and the curvature tensor K a p is 

K a p = N • D a D IL = R a K aP 

1 ( {l + p)-dlp + 2{de P f 



K, 



a/3 



y/1 + (Vp) 2 V - sin 9dg §fi + 2d e pd 4> p (1 + p) sin 2 9 - sin 9 cos 9d e p - d 2 p + 2(9 p) i 



- S m6d e ^+2dgpd^p 



(5.64) 
(5.65) 

(5.66) 



(5.67) 
(5.68) 



where (Vp) 2 = (dgp) 2 + „ (d^p) 2 . The Hamiltonian is a functional of and p(f2) at the fixed area A. To 

find the equilibrium form of ip(£l) and p(f2), we need to minimize H over ip(Q) and p(O). There is a considerable 
simplification if we restrict our attention to the neighborhood of the transition temperature Tq. Near Tq, to order ip 4 , 
the variation of the curvature energy is of order up 2 . When this is comparable to the Landau-Ginzburg free energy 
for ip, we only need to keep couplings of order pip 2 . To this order, the Landau-Ginzburg Hamiltonian for ip becomes 



i M |V| 4 ) + X -K n J d£lg a P(d a - inA° a )iP(d + mA%)r 



+ J dttlp (r\iP\ 2 + iu|V| 4 ^) + J dnj a 5A a , 



(5.69) 



where A° e = and A , 



■ cos9. In this formula, 



-1/2 



SH 

SA a 



^f^(^d r - pdf,ij> + 2mA°|V>| 2 ) 



(5.70) 



has been evaluated at p; m = 0, Ro = R and 5A a — j^dpp. In terms of p, the curvature energy can be written as 



n K = l -n j dn((w 2 + 2) P y 



where 



V 2 p 



1 



Hence, to order ip 4 , we can write H as H = H + 5H with 



sinf 

- -wo 



-dg sin 9dgp 



1 



sin 2 9 



(5.71) 



(5.72) 



dn[r\^\ 2 + -u\^\ i 



-I 

X -K n J dnf a f\d a - inA° a )iP(dp + inA°)r, 



(5.73) 
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which does not depend on p(O) and 

+ i K y^((V 2 + 2)p) 2 , (5.74) 
which contains p(SY). This can be rewritten as 

h = n Q + J dnp(fi)$(n) + \* j dmrt p {Q)m{q,, n') P (n') (5.75) 



with 



= ^ (r#| 2 + \u\^ - -I_d Q ( 7 ^) (5.76) 



OO I 



m(q, n') = £ H y ; m *( fi )(J - 1) 2 (; + 2) 2 r ; m (r!'). (5.77) 



(=2 m=- 

Minimizing over p(ft), we obtain 



P (n) = -I / dfi'M-^n.n')*^'). ( 5 - 78 ) 

where M -1 (f2, f2') satisfies 

y rf0 3 M(0!, n 3 )M- 1 (0 3 , fi 2 ) - <5(i>2)(fii - ^2) (5.79) 

and 5(;>2)(^i — ^2) is the Dirac delta function in spherical harmonics space with I > 2. By substituting this for p(Q) 
into 7i, we have 

+ ±K n J dSlfv^da - inA° a )iP(dp + inA%)p 

1 00 ^ | 5L |2 



where 



= J dOF; m *(0)$(^) (5.81) 
,4r 1 

$(0) = ^H 2 -— 9 Q ( 7 ^). (5.82) 

Before proceeding with our analysis of the n-atic transition on a deformable sphere, it is useful to recall Abrikosov's 
calculation of the transition to the vortex state. The first step is to calculate the eigenf unctions of the harmonic 
part of Hgl when V x A = H. These can be divided into highly degenerate sets separated by an energy gap, 
huj c = 2Ce*H. In the Landau gauge with A = (0, Hx,0), the eigenfunctions in the lowest energy manifold are 
ijjk = e lhy e~ e H ( x ~ x h) w here Xk = k/e*H. The order parameter ip(x) of the ordered state is expressed as a linear 
combination tp{x) = ^ CfcV'fc; where the complex parameters Ck are determined by minimization of 7Ygl- 

With this analogy in mind, we will follow Abrikosov's treatment of the superconducting transition near Hq 2 to 
study the development of n-atic order on a sphere. We first diagonalize H in the harmonic level when p — and 



13 



Ag = Ag, A^ — that is, we determine the functions ip that satisfy KAD a D a ip = eip, for p = 0, A Q = 
Nonlinearities arise from the | ip | 4 term in Landau- Ginzburg Hamiltonian for ip and from the fact that A a and 
g a (i depend nonlinearly on the deviation from the ideal spherical shape. To minimize H c s, first we seek the lowest 
energy configurations of the operators corresponding to the harmonic terms only. Then, we linearly combine these 
eigenfunctions to get the function which has the lowest energy for 7Y c ff • That is, the manifold of the lowest energy 
eigenstates of the harmonic terms is highly degenerate. Nonlinear terms pick out the lowest energy state which is the 
linear combination of the degenerate eigenstates of the harmonic terms and resolve the degeneracy. We divide Tt e g 
into the harmonic part, 7iharj and the nonlinear part, H n \: 



= ^ / dnr \^\ 2 + \ R n J d ^f aP \d a - inA° a )iP(dp + ^nA° )V* 



The differential operator for Hhai is 



where 



=2 m=-l 



(Z - 1) 2 (Z + 2)2 



1 SHhar , A K n 



(5.83) 
(5.84) 

(5.85) 



A = - 



(d a -inAl)^ff a(3 (dp-inAl). 



In the spherical polar representation, A has the form 

1 



A 



sin 6 



« i ■ o« \ , 1 «2 2 cos 2 6» cos# 
smfl v sinO sin 6* 



One possible form of the eigenfunctions of A has the structure of sin™ 9. 

A sin™ 9 = nsm n 6. 



(5.86) 



(5.87) 



(5.88) 



The complete spectrum and the eigenfunctions are derived in the Appendix. Introducing the projection representation, 
we can reparametrize the sphere using a (z, z) coordinate system defined by 



9 

z = tan -e 
2 



i(f> 



z = tan -e 
2 



-i(f> 



The inverse transformation gives 



9 = 2 tan 1 \fzz 



2i z 



(5.89) 



(5.90) 



The corresponding metric tensor, spin connection, and the differential operator, A, in this representation can be 
written as follows 



1 

{l + \z\ 2 V \ 1 



(5.91) 



i i 



M 2 



2iz 1 



1 1 



M 5 



2iz 1 



(5.92) 



A = -(1 + \z\ 2 f[{d z d- z - A° X A%) - 2in(A% + A%)}. 
The eigenfunction sin™ 9 can be written as 

2\z\ 



f(\*n = 



(5.93) 



(5.94) 
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This eigenfunction has vortices at \z\ = and \z\ = oo. Now, the vortices can be arbitrarily placed by multiplying 
Q(z) which is a function of z alone. From the differential operator A, we see 

A/(|z| 2 ) = n/(|z| 2 ) 
A/(|z| 2 )Q(z) = n/(|z| 2 )Q(z). (5.95) 

Thus, the full eigenfunctions of A are the product of /(|z| 2 ) and the arbitrary function Q(z). In particular, for n = 1, 
we have 

/(M 2 ) = (5-96) 



By taking Q(z) 



(z-z 1 )(z-z 2 ) 



i in <-\ 2 I Z I {z-z l ){z-z 2 ) , 

ip = Const.) , 5.97) 

1 + \zr z 



which means that ip has the vortices at z = Z\ and z — z 2 and that ip goes to z 1 z 2 e -4 ^ as z goes 0, ip goes to e 1 ^ as 
z goes to oo. This means there is an enormous degeneracy in the lowest energy manifold of the harmonic terms. In 
the lowest energy manifold of the harmonic terms, tp will have exactly 2n zeros specifying the vortex positions, each 
of strength 1/n. 

In view of this vortex strength and distribution, we choose Q{z) as 



2n 

a 



Q(*) = ^rII(*-*), (5-98) 

i=l 

where a, Zi are arbitrary and Zj ^ Zj if i 7^ J. By minimizing 7Y c ff with the degenerate eigenfunctions of A, 

/ 911 \ n 1 ?" 

^ = a (Trf^) -11^-^, (5-99) 



the non-linear terms pick out the set {zi | 1 < i < 2n} that gives the lowest energy for TC c s ■ This set {zi} also gives 
\ip\, vortex positions, and the shape change R of the membrane. In terms of spherical coordinates, 



ip = (Const.) [J sin - cos _L e *(^-^)/2 _ cos _ sin _* e -i«>-^)/2 _ ^ p(O) (5.100) 
i=1 V ^ 2 2 2 / 

and 

/ 2 2n 

l^| 2 -T7ll( 1 - C0S ^) ( 5 - 101 ) 



2 — 1 

where cos7i = cos6>cos6>i + sin 8 sin 0i cos(0 — (pi), = <Pi) specifies the position of the i-th zero of ip, and AT is 
the normalization factor which makes 



dn\ip\ 2 = ipl. (5.102) 

since ip of Eq. (|5.100D 

is the most general function in the lowest energy manifold, the order parameter and vesicle 
shape just below the transition temperature are obtained by minimizing 7Y e ff over ipg and the positions of zeros. 
Inserting this expression for ip, the free energy becomes 



dil ({r-r c )\iP\ 2 + ^u\iP\^j 



1 00 ^ I rf\ 1 2 
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where r c = —AtthKa/ A and 



= -nif I ^ 2 J|(l-cos 7j ) 



1 + Sn 



cos 74 



r^(l-cos7i) ^ (1- cos 7^(1 -cos 7,) 



cos "fij — cos 7i cos t,' 



Minimization over ?/>o leads to -00 = — K{^i}]/u[{f2i}] ancl the effective free energy density 

/[i " i|J 2«[{fi i }]' 
depending only on the positions {f^} of the zeros of i^iuf) where 



■[{Hi}] = (r-r c ) / dft|P(ft)| 2 
7 



oo / 
Z=2 m=-i v 



|4> 



lm | 



i) 2 a + 2) 2 



(5.104) 
(5.105) 

(5.106) 
(5.107) 



The next step is to minimize /[{Oj}] over £1^. Shape changes are described by p(Q) = Et2E m =-ift'"'m(^) 
where 



1 



Pin 



k{1-\) 2 {1 + 2) 2 
k{1- 1) 2 (Z + 2) 2 



dQY; n *{n)^{VL) 



*j m ({fii}), 



(5.108) 



with {Qi} that minimize In the disordered phase above the transition temperature, ipo = 0. Hence, p = 

and i?o = R- The vesicle shape is spherical. In the ordered phase below the transition temperature, we have been 
able to evaluate /[{fii}] analytically for n — 1 and n — 2. For these two cases, we find as expected that the minimum 
energy configurations are, respectively, those with zeros of ip(Sl) at opposite poles and at the vertices of a tetrahedron. 
The shape function p(Q) = {nK n / 'K)tpQp^ n ' (f2) associated with n-atic order is proportional to tpQ ~ (r — r c ) to the 
order of our calculations. In general, the Legendre decomposition of p( n ' will contain Legendre polynomials of order 
2n. For n = 1 and n = 2, we find 



p«(0) = 



(5.109) 



H y ' 135 V 7 



Y 3 °(n)-\/-(Yi(n)-Yf 3 (n)) 



3645 



(n 3 (^) 



(5.110) 



Figure [l] shows the shapes described by these functions. The transformations from the initial spherical shape to 
distorted shapes occur continuously. Our calculations for the shape are valid to order (r — r c ). As temperature is 
lowered, higher-order terms in (r — r c ) and higher-order spherical harmonics are needed to describe the equilibrium 
shape. In Ref . , a variational function for ip and spherical harmonics up to order 8 were used to calculate the shape 
for n = 1 for temperatures well below the transition. We use symmetry considerations to treat n = 3,4, and 6. For 
n = 3 and n = 6, we expect the zeros of to lie, respectively, at the vertices of an octahedron and an icosahedron. 
We find the shape functions for n = 3 and n = 6 



Y 4 °(n) + J-(Yi(n) + Y 4 ~ 4 (n)) 



60060 



'137T 



F 6 °(0) 



-(y 6 4 (tt) + y 6 - 



(5.111) 
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15488 r- 
275559375 ' 

256 
+ 692803125 

512 

"6834953125 1 



Y 6 °(n) + x l-(Yi(n)-Y e - 5 (n)) 



1 R7 

^(ft) - Y^ft)) + -J-(Y™(Sl) + Yr 10 (n)) 



(5.112) 



Figure shows the shapes described by these functions. For n = 4, following the calculations of Palffy-Muhoray 
plfl , we expect the zeros of tp(Q) to lie at the vertices of a distorted cube obtained by rotating its top face about its 
fourfold axis by 7r/4 and compressing opposite faces. For this case, we minimized the energy over a single parameter 
determining the separation between opposite rotated faces and obtain 



pW(fi) = -0.0184631F 2 °(ft) 
+0.0201099 Y 6 ° (ft) 
-0.00159671 [y 8 ° - 



- 0.305216F 4 (ft)0.0335485(F 5 4 (ft) 

- 0.00721177(r7 t (ft) + Yf 4 (ft)) 
0.578489(F 8 8 (ft) +F t f 8 (ft))] . 



Y 5 - 4 m 



(5.113) 



Figure || shows the shape for n = 4. 



VI. CONCLUSIONS 



We have presented the general n-atic Hamiltonian in terms of the n-th rank spherical tensors and have shown 
there exists an effective KT transition for n 2 K n / 'k C 1/4. Also we have presented here an analysis of the mean-field 
transition to n-atic order on a fixed-area surface of genus zero and the corresponding continuous shape changes below 
the KT transition temperature in this parameter region. We are really dealing with two kinds of order: n-atic order 
and the positional order of vortices. In mean-field theory, these two kinds of order develop simultaneously. Our 
analysis is very similar to that of Abrikosov for the transition from a normal metal to a vortex lattice in a type II 
superconductor at Hq 2 and ignores the effects of fluctuations. In superconductors, fluctuations drive the normal-to- 
superconducting transition first order in a nonvanishing magnetic field |22|. In the Abrikosov phase, fluctuations of the 
vortex lattice destroy superconductivity but not long-range periodic order. In two-dimensions, screening of vortices 
drives the Kosterlitz-Thouless transition in an infinite superconductor in zero field to zero. Both of the above effects 
may be important for n-atic order on a sphere. We have been able to find in the parameter region n 2 K n / n <C 1/4, the 
effect of shape fluctuations is irrelavent and will not lead to qualitative changes in our results. However, outside this 
region, i.e. for n 2 K n /K 3> 1/4, shape fluctuations are relavent. The interaction between the vortices is screened to 
become the massive sine-Gordon theory and hence is equivalent to a neutral Yukawa gas on a sphere. Consequently, 
this screening effect derives the Kosterlitz-Thouless transition temperature to zero. The vortices are always unbound 
for a non-zero temperature. 

In Ref. (2^] , Evans has discussed the shape changes of the n-atic Hamiltonian in the opposite limit (high temperature 
limit) to the mean-field limit described in this paper. Using the lowest Landau level approximation, the effect 
of thermal fluctuations are discussed. Although we agree with him in that the amplitudes of deformations have 
zero thermal average, since the lowest Landau levels include only the ground state vortices we are not sure if this 
approximation is valid in the high temperature phase where thermally excited vortex-antivortex pairs are always 
unbound. Thus we believe the effects of fluctuations to the shape changes in the disordered state deserve further 
investigation. 



APPENDIX A: COMPLETE SPECTRUM OF A 



We present the complete spectrum of the differential operator, A, for the harmonic part, Hhar- The operator A in 
the spherical polar representation, 



A = — ^— (desmOdg + -t^(<9<a + incos9) 2 ) , 
smv \ smO J 

is similar to the differential operator —(V — iA) 2 for the superconductivity where 



(1.1) 
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V = 9d 



smt 



(1.2) 



and 



A = 9Ao 



(1.3) 



The Schrodinger equation on a sphere in the field of a Dirac magnetic monopole of strength g at the origin plus an 
infinitely long and thin solenoid carrying flux F along the z axis is described by 



1 



(V - iA) 2 il> = Eip 



with 



A e = 0; = _£_(i - cosfl) + — — 
sm 9 2ir sm ( 



(1.4) 



(1.5) 



The monopole has been chosen to have a Dirac string singularity along the negative z axis. Setting g = n, F = —2irn, 
we find 



A e = ; A 4 



ncosQ 



suit 



(1.6) 



and — (V — iA) 2 is exactly identical to A. Hence, the differential operator, A, is equivalent to the Hamiltonian 
describing the field of a Dirac magnetic monopole of strength n at the origin plus an infinitely long and thin solenoid 
carrying flux 1-kti along the positive z axis. Following the method in Ref. |24], we define the formal differential operator 



J = — ir x (V — iA) — nr 



with components 



,±id 



±dg + i cot 9{d ( j > 



J± = J x ± iJ y = e 
J z = —i{d<f, + in) - n 
which obeys the formal commutation relations 

[J,W] = ; [Ji, Jj] — ieijkJk-, 



nsm9 



1 + cos ( 



for sin 6 ^ 0. We find 



Let us denote 



where we set 



Then 



J 2 = n 2 + A 
=> A = J 2 -n 2 . 



j 2 *(^) =i(i + i)*(n), 



= exp[im(j)}P (cos 9) ; m = 0, ±1, ±2, • 



The eigenvalues of J 2 are found to be given by 1(1 + 1) with 

I = k+ i (\m + n\ + \m - n\) ; k = 0, 1, 2,3, • • • ; m = 0, ±1, ±2, • 



(1.7) 

(1.8) 
(1.9) 

(1.10) 

(1.11) 
(1.12) 
(1.13) 
(1.14) 

(1.15) 



The corresponding orthonormal cigenfunctions are 
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^(n) = C km ' 2 ' ( 1± |^) ' 2 ' • Pi |m+ " l>l -"'W)e™^ (1.16) 

where p^ a ' b ^ denote Jacobi polynomials, and 

_ / fc!r(|m + n\ + \m-n\ + k + l)T(\m + n\ + |m - n\ + 2fc + 1) \ 5 
fem_ ^ 47rr(|m + n| + fc + l)r(|m - n| + fc + 1) / ' ( ' j 

The complete spectrum, <5, of the operator A is 

<5 = (fc + — (|to + n| + \m — n\))(k + 1 + — (\m + n\ + |m — n|)) — n 2 

fc = 0,l,2,3,--- ; m = 0,±l,±2,.--. (1.18) 

5=(n + k) 2 -n 2 + (n + k). (1.19) 

<5 = (|to| + fc) 2 -n 2 + (|to| + fc). (1.20) 



For |to| < n, 



For ItoI > n, 



Thus, 

5 = (max(n, \m |) + fc)(max(n, |to|) + fc + 1) - ?i 2 . (1.21) 
The ground state energy is generated when fc = 0, n > |to|. 

5 = n(n + l)-n 2 =n (1.22) 
The corresponding ground state orthonormal eigcnfunction has the form 

¥°' m \n) = C 0m (sm-) (cos^j e™* (1.23) 



which is the same as that in Eq. ( |5.100Q 
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TABLE I. The nth rank symmetric traceless tensors, Q' n ', for n = 1, 2, 3, 4, and 6 in 2-dimensional space. 



— Trial rn a 2 ^l& 3 l^(,8aia 2 iri a3 8a 2 a 3 7Tl a ± ^ a 3 a l ) 



}(3) 

; a.i Q2 C13 



^(8 aia2 m a3 m a4 + 5 more terms by permutations) 
+ j4(5a 1 a 2 fia 3 a i + 2 more terms by permutations) 
Qa}a 2 a 3 a 4 a 5 a 6 = m ai m a2 m a3 m a4 m as m a6 - 1 \(8 aia2 m a3 m a4 m a5 m a6 + 14 more terms by permutations) 



+ -^(8 aia2 S a3a4 m a5 m a6 +44 more terms by permutations) 
77 (8 ai g 2 8 a3Cli 8 a5a6 + 14 more terms by permutations) 



m 



Q a 

Tr^ai r ffl j ct^ Qa.i a 2 »3 4 (^a 1 Q2 Q03 H~ ^a2 Q 3 Qcx\ H - ^0301 Qa.2 ) 

m ai m a2 m a3 m a4 = <5aia 2 a 3 a 4 + gOSaic^Qas^ + 5 more terms by permutations) 

+ |(5 c « 1 a2 | 5a3c«4 +2 more terms by permutations) 
m ai m a2 m a3 m a4 rn a5 rn a6 = Qai a203(I4<15a6 + ^((Ja^Q^^sae + 14 more terms by permutations) 

4-^(i5a 1 a2 8a 3 a 4 Qa2a 6 +44 more terms by permutations) 
+ jg{8aia 2 8a 3 a i 8 a5 a 6 + 14 more terms by permutations) 
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FIG. 1. Mean-field shapes of deformable surfaces of genus zero with vector order (n = 1) on the left and nematic order 
(n = 2) on the right. Above the mean-field transition temperature, the equilibrium shape is spherical for all n. Below the 
transition temperature, the equilibrium shape depends on n and has a ellipsoidal form with 2 ground state vortices with strength 
1 located at opposite poles for n = 1 and a tetrahedral form with 4 ground state vortices with strength 1/2 located at the 
vertices of a tetrahedron for n — 2. 



FIG. 2. Mean-field shapes of deformable surfaces of genus zero with tertic order (n = 3) on the left and hexatic order (n = 6) 
on the right. Below the transition temperature, the equilibrium shape has a octahedral form with 6 ground state vortices 
with strength 1/3 located at the vertices of a octahedron for n — 3 and a icosahedral form with 12 ground state vortices with 
strength 1/6 located at the vertices of a icosahedron for n = 6. 
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FIG. 3. Mean- field shapes of deformable surfaces of genus zero with quartic order (n — 4). Below the transition temperature, 
the equilibrium shape has a distorted cubic form with 8 ground state vortices with strength 1/4 located at the vetices of a 
distorted cube obtained by rotating its top face about its 4-fold axis by 7r/4 and compressing opposite faces. 
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